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ABSTRACT 

It is common to estimate black hole abundances by using a measured correlation between black hole 
mass and another more easily measured observable such as the velocity dispersion or luminosity of 
the surrounding bulge. The correlation is used to transform the distribution of the observable into 
an estimate of the distribution of black hole masses. However, different observables provide different 
estimates: the M, — a relation predicts fewer massive black holes than does the M, — L relation. 
This is because the a — L relation in black hole samples currently available is inconsistent with that 
in the SDSS sample, from which the distributions of L or tr are based: the black hole samples have 
smaller L for a given a or have larger a for a given L. This is true whether L is estimated in the 
optical or in the NIR. The a~ L relation in the black hole samples is similarly discrepant with that in 
the ENEAR sample of nearby early-type galaxies. This suggests that current black hole samples are 
biased towards objects with abnormally large velocity dispersions for their luminosities. If this is a 
selection rather than physical effect, then the Af, — a and Af, — L relations currently in the literature 
are also biased from their true values. We provide a framework for describing the effect of this bias. 
We then combine it with a model of the bias to make an estimate of the true intrinsic relations. While 
we do not claim to have understood the source of the bias, our simple model is able to reproduce 
the observed trends. If we have correctly modeled the selection effect, then our analysis suggests that 
the bias in the (Af, |tT) relation is likely to be small, whereas the (Af, |i) relation is biased towards 
predicting more massive black holes for a given luminosity. In addition, it is likely that the Af, — L 
relation is entirely a consequence of more fundamental relations between Af, and a, and between a 
and L. The intrinsic relation we find suggests that at fixed luminosity, older galaxies tend to host 
more massive black holes. 

Subject headings: galaxies: elliptical — galaxies: fundamental parameters — black hole physics 



1. INTRODUCTION 

Several groups have noted that galaxy formation and 
supermassive black hole growth should be related, and 
many have modeled the joint evolution of quasars and 
galaxies (e.g., Monaco et al. 2000; Kauffmann & 
Haehnelt 2001; Granato et al. 2001; Cavaliere & Vit- 
torini 2002; Cattaneo & Bernardi 2003; Haiman et al. 
2004; Hopkins et al. 2004; Lapi et al. 2006; Haiman et 
al. 2006 and references therein). However, the number of 
black hole detections to date is less than fifty; this small 
number prevents a direct estimate of the black hole mass 
function. On the other hand, correlations between Af, 
and other observables can be measured quite reliably, if 
one assumes that the Af,— observable relation is a sin- 
gle power-law. Af, is observed to correlate strongly and 
tightly with the velocity dispersion of the surrounding 
bulge (Ferrarese & Merritt 2000; Gebhardt et al. 2000; 
Tremaine et al. 2002), as well as with bulge luminosity 
(McLure & Dunlop 2002) and bulge stellar mass (Mar- 
coni & Hunt 2003; Haring & Rix 2004): in all cases, a 
single power law appears to be an adequate description. 

Since it is considerably easier to measure bulge prop- 
erties than Af, , black hole abundances are currently esti- 
mated by transforming the observed abundances of such 
secondary indicators using the observed correlations with 
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Af, (Yu & Tremaine 2002; AUer & Richstone 2002; Mar- 
coni et al. 2004; Shankar et al. 2004; McLure & Dunlop 
2004). Unfortunately, combining the Af,— cr relation with 
the observed distribution of velocity dispersions results 
in an estimate of the number density of black holes with 
Af, > lO^Af© which is considerably smaller than the es- 
timate based on combining the Af, — L relation with the 
observed luminosity function (Tundo et al. 2006). 

The fundamental cause of the discrepancy is that the 
luminosity and velocity dispersion functions are obtained 
from the SDSS (Blanton et al. 2003; Sheth et al. 2003), 
but the a—L relation in the SDSS (Bernardi et al. 2003b) 
is different from that in the black hole samples (e.g. Yu 
& Tremaine 2002). In Section^ we argue that the SDSS 
a — L relation is not to blame, suggesting that it is the 
black hole samples which are biased. 

If black hole samples are indeed biased, then this bias 
may be physical or it may simply be a selection effect. 
If physical, then only a special set of galaxies host black 
holes, and the entire approach above (of transforming 
the luminosity or velocity dispersion function) is compro- 
mised. On the other hand, if it is a selection effect, then 
one is led to ask if the Af,— observable relations currently 
in the literature are biased by this selection or not. The 
main purpose of this work is to provide a framework for 
describing these biases. Section [3] describes Monte-Carlo 
simulations of the selection effect. It shows an example 
of the bias which results, and provides an analytic model 
of the effect. Our findings arc summarized in Section |4l 
where estimates of the intrinsic black hole mass function 
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Fig. 1. — Comparison of the a — L relation in the black hole sample of Haring & Rix (2004) (symbols) with that from local early- type 
galaxy samples (lines). We provide two estimates of the bulge luminosity for the Haring & Rix compilation: the open circles come from 
rescaling the bulge luminosity to the SDSS r-band, and the filled circles from assuming a mass-to-light ratio to scale the bulge masses 
Mbuigg to luminosities. Dotted line shows a fit to the filled circles. Solid and short dashed lines show {<t\L) in the SDSS-B06 and ENEAR 
samples. Dot-dashed line shows the biased relation obtained from the ENEAR sample if one fails to account for the fact that the velocity 
dispersions played an important role in determining the distances from which the luminosities were estimated. The black hole compilation 
is clearly offset from the SDSS-B06, the ENEAR and even the biased ENEAR relations. Appendix A gives the fits shown, and describes 
exactly how the black hole sample was compiled. 



4>{M,), and the intrinsic M, — a and M, — L relations 

are provided. 

Throughout, we assume a spatially flat background 
cosmology with = 0.3 and a cosmological constant. 
All luminosities and masses have been scaled to Hq = 
70 km s~^ Mpc~^ 

2. MOTIVATION: THE a-L RELATION 

This section shows that the a — L relation of early- 
type galaxies is significantly different from that in current 
black hole samples. This difference exists whether L is 
computed in a visual band, or in the near-infrared. 



2.1. Remarks on SDSS photometry and spectroscopy 

Since 2002, when the Bernardi et al. (2003a) sam- 
ple was compiled, a number of systematic problems with 

SDSS photometry and spectroscopy have been discov- 
ered. The initial (SDSS-DRl and earlier) photometry 
overestimated the luminosities and sizes of extended ob- 
jects. Although later releases correct for these problems, 
they still overestimate the sky level in crowded fields 
(e.g. Mandelbaum et al. 2005; Bernardi et al. 2006; 
Lauer et al. 2006), and so they tend to underestimate 
the luminosities of extended objects in such fields. This 
problem is particularly severe for early-type galaxies in 
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nearby clusters. In addition, the SDSS slightly overes- 
timates the velocity dispersions at small a (the velocity 
dispersions in Bernardi et al. 2003a differ from the val- 
ues reported in the SDSS database and are not biased 
at small a). Bernardi (2006) discusses an analysis of the 
a — L which results once all these systematic effects have 
been corrected- for: the net effect does not significantly 
change the a — L relation reported by Bernardi et al. 
(2GG3b). Nevertheless, the results which follow are based 
on luminosities and velocity dispersions which have been 
corrected for all these effects. Because they are not eas- 
ily obtained from the SDSS database, in what follows, 
we refer to the corrected SDSS sample as the SDSS-B06 
sample. 

Bernardi (2006) shows that the a — L relation in the 
SDSS-B06 sample is consistent with that in ENEAR, the 
definitive sample of nearby (within 7000 km s^^) early- 
type galaxies (da Costa et al. 2000; Bernardi et al. 2002; 
Alonso et al. 2003; Wegner et al. 2003), provided that 
one accounts for the bias which arises from the fact that 
fj played an important role in determining L. Hence, the 
local (ENEAR) and not so local determinations (SDSS- 
B06) of the (7—L relation are consistent with one another. 

2.2. The a — L relations in the optical 

Figure [T] compares the correlation between a and L in 
ENEAR and SDSS-B06 with that in the recent compila- 
tions of black holes by Haring & Rix (2004) . The dotted 
line shows a fit to {cr\L) defined by the filled circles. Ap- 
pendix A describes this and other black hole compilations 
we have consider in the following figures. It also gives the 
fits shown in this and the other figures. 

The solid line shows {cr\L) in the SDSS-B06 sample. 
The dot-dashed line shows a determination of the EN- 
EAR a — L relation which attempts to correct for pe- 
culiar velocity effects in the distance estimate. This 
correction is known to bias the (cr|L) relation towards 
a steeper slope (Bernardi 2006). The short-dashed line 
shows the result of using the redshift as distance indica- 
tor when computing the ENEAR luminosities; while not 
ideal, comparison of this relation with the dot-dashed 
line provides some indication of the magnitude of the 
bias. Note in particular that using the redshift as dis- 
tance indicator brings the ENEAR relation substantially 
closer to the SDSS-B06 relation, for which the redshift is 
an excellent distance indicator. 

The agreement between ENEAR and SDSS-B06 sug- 
gests that the difference between the a — L relation in 
the black hole sample (dotted) and in SDSS-B06 (sohd) 
cannot be attributed to problems with the SDSS deter- 
mination. 

The agreement between ENEAR and SDSS-B06, and 
the disagreement with the black hole samples strongly 
suggests that the a — L relation in black hole samples is 
biased to larger a for given L, or to smaller L for a given 
(7. Similar conclusions are obtained if we use Kormendy 
& Gebhardt (2001) or Ferrarese & Ford (2005) compi- 
lations. At a given luminosity, the black hole samples 
have logcr larger than the ENEAR or SDSS-B06 rela- 
tions by ~ 0.1 dex. As a result of this difference, smaller 
luminosities can apparently produce black holes of larger 
mass than do their associated velocity dispersions (e.g., 
Yu & Tremaine 2002). Tundo et al. (2006) show this in 
some detail. 



2.3. The a ~ L relations in the near-infrared 

We have considered the possibility that the discrep- 
ancy between the a — L relations may depend on wave- 
length. Figure [2] presents a similar analysis, but now the 
luminosities are from the i^-band. In the case of the 
early-type sample, the SDSS-B06 sample was matched 
to 2MASS, and then the 2MASS isophotal magnitudes 
were brightened by 0.2 mags to make them more like to- 
tal magnitudes (Kochanek et al. 2001; Marconi & Hunt 
2003). The solid line shows the resulting (ctIXk) relation. 
It is very similar to that in the r— band ii r — K = 2.8. 

Marconi & Hunt (2003) matched their black hole com- 
pilation to 2MASS, and then recomputed the photometry 
to determine Lk- The open triangles show the a — Lk 
relation for the subset of objects which are in the Haring 
& Rix sample. The filled circles are obtained by tak- 
ing these objects (i.e. those in both Marconi & Hunt 
and Haring & Rix), and then shifting the Mbuigc based 
r— band luminosities to K assuming r — K — 2.8. Dashed 
and dotted lines show fits to these relations. Compari- 
son with the solid line (2MASS-B06) shows a significant 
offset — the black hole samples are different from the pop- 
ulation of normal early-type galaxies. 

3. BIASES IN THE Af. - a AND ill. - L RELATIONS 

If black hole samples are indeed biased, then this bias 
may be physical or simply a selection effect. One is then 
led to ask if the relations currently in the literature are 
biased by this selection or not. To address this, we must 
be able to quantify how much of the M, — a correlation 
arises from the M, — L and L^a relations, and similarly 
for the M, — L relation. 

3.1. Demonstration using mock catalogs 

Our approach is to make a model for the true joint 
distribution of M,, L and cr, and to study the effects 
of selection by measuring the various pairwise correla- 
tions before and after applying the selection procedure 
to this joint distribution. In what follows, we will use 
V = log(cr/km s~^), and with some abuse of notation, 
M, to denote the logarithm of the black hole mass when 
expressed in units of Mq and L to denote the absolute 
magnitude (since, in the current context, M for absolute 
magnitude is too easily confused with mass). We assume 
that correlations between the logarithms of these quanti- 
ties define simple linear relations, so our model has nine 
free parameters: these may be thought of as the slope, 
zero-point and scatter of each of the three pairwise cor- 
relations, or as the means and variances of L, V and M, 
and the three correlation coefficients rvL, f»v and r.^. 

To illustrate, let (V) and cry denote the mean and rms 
values of V, and use similar notation for the other two 
variables. Then the mean y as a function of L satisfies 

(Tl 

and the rms scatter around this relation is 



crv\L ^ crv yl- rl^L- (2) 

The combination ay|^ = rvL <^v/<^l is sometimes called 
the slope of the {V\L) relation, and {V) — ay^i^ (L) is the 
zero-point. Similarly, the slopes of {M,\V) and (A/, |L) 
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Fig. 2. — Same as Figure 1 but now using the X-band luminosity. Triangles show the subset of the Haring & Rix sample for which X-band 
luminosity estimates are available (from Marconi & Hunt 2003). Solid circles show the result of rescaling the Afbuigp-based luminosities of 
these objects to the ii'-band using r — K = 2.8. Dotted and dashed lines show fits to the ((t|Lk) relations defined by these points. Solid 
line shows ((t|Lk) obtained by matching the SDSS-B06 early-type sample to the 2MASS database, and then using the ii'-band luminosities 
with the SDSS-B06 velocity dispersions. This fit is very similar to the solid line shown in the previous panels '\l r — K = 2.8. These fits are 
given in Appendix A. 



are a,|y = r.ycr./cry and a,|i = r,Lcr,/aL respec- 
tively. Notice that the zero points are trivial if one works 
with quantities from which the mean value has been sub- 
tracted. We will do so unless otherwise specified. 

Our analysis is simplified by the fact that five of the 
nine free parameters are known: the mean and rms values 
of L and V as well as the correlation between V and 
L have been derived from the SDSS-B06 sample (these 
values are similar to those reported by Bernardi et al. 
2003b, who also show that the distribution of V given L 
is Gaussian). Hence, we use 

CTL = 0.84, cry = 0.11, and tvl = -0.78, 

to make a mock catalog of the joint V — L distribution 
in the SDSS-B06 early-type galaxy sample. Sheth et 



al. (2003) show explicitly that the resulting mock cata- 
log provides luminosity and velocity dispersion functions 
which are consistent with those seen in the SDSS. 

We turn this into a catalog of black hole masses by 
assuming that, like V at fixed i, the distribution of M, 
at fixed V and L is Gaussian. Hence, for each V and L, 
we generate a value of Af, that is drawn from a Gaussian 
distribution with mean 

{M,\L,V) _ (M.) ^ L - (L) r,L - r^yryL 
a, a, (TL '^-TvL 

V - {V) r,v - r.LrvL , , 

ay 1 - ri^r 
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Fig. 3. — Comparison of the joint distribution of M,, a and L, before (open symbols) and after (filled symbols) applying a selection cut. 
Dashed and solid lines show the intrinsic and observed correlations. The probability of selection is assumed to increase as one moves to 
the right of the dashed line in the bottom left panel. 



and rms 
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(4) 



•VL 



(see, e.g., Bernardi et al. 2003b). The combination 
"i'tvryi^ measures how much of the correlation between 
M, and L is simply a consequence of the correlations 
between M, and V and between V and L. The combi- 
nation TtLryL can be interpreted similarly. Notice again 
that the mean values of the three quantities mainly serve 
to set the zero-points of the correlations. Since we are 
working with quantities from which the mean has been 
subtracted, our model for the true relations is fixed by 
specifying the values of a,, r,y and r.^. 

To study how selection effects might have biased the 
observed M, — u and M, — L relations we generate a mock 



galaxy -l-black hole catalog as described above, model the 
selection effect, and then measure the three pairwise cor- 
relations in the selection-biased catalog. Figure [3] shows 
an example of the result when 

a, = 0.49, r,v = 0.88, and r.L = -0.70. 

We have tried a variety of ways to specify a selec- 
tion effect. In Figure [3l we supposed that no objects 
with luminosities brighter than {L\V) were selected, and 
that the fraction of fainter objects which are selected is 
erf(a;/V2), where x = {L - {L\V))/aL\v. The open di- 
amonds in each panel show the objects in the original 
catalog, and the filled circles show the objects which sur- 
vived the selection procedure. The panels also show the 
slopes of the relations before (dashed) and after (solid) 
selection. Note that because the selection was made at 
fixed V, the slope of (M, |y) in the selected sample is 



6 



A 



C71 

o 



V 
I 

I 
m 



1.0 



0.5 ■ 



0.0 



-0.5 - 



-1.0 




-0.2-0.1 0.0 0.1 0.2 
Log,o <T-<Log,o o\Mf> 



1.0 



0.5 



I 

m 



V 0.0 
I 

r 

OD 

^ -0.5 
-1.0 




-0.2-0.1 0.0 0.1 0.2 
Log,o CT-<Logio o\M,> 



V 
I 



1.0 



0.5 - 



0.0 



-0.5 ■ 



-1.0 



1 1 

o 
• 

o 

• 

• 

O ^ 


o • 

o 

o 


o. go* \iS 

o S « 


» o f 

1 1 



-0.2-0.1 0.0 0.1 0.2 
Log,o a - <Log,o o> 



1.0 



0.5 



V 0.0 
I 

r 

OQ 

^ -0.5 



-1.0 



1 

o 

<i 

o 


1 

• o 

o •« 

o» « , o 
^ $o , oo » 




■ 


F.-^ o 

Q 

o 

I 



-2-10 1 
M,-<MjLog,o o> 



Fig. 4. — Residuals from the A/, — cr relation do not correlate with residuals from the a — L correlation (top left); whereas residuals 
from the Af, — L relation do (top right). They also correlate with residuals from the L — cr relation (bottom right), and weakly with a 
(bottom left). In all pan els, open circles show the intrinsic correlations and filled circles show correlations after applying the selection cut. 
Comparison with Figure lA2l shows that the filled circles trace similar locii to those seen in the Haring & Rix (2004) compilation. 



not changed from its original value. However, {M,\L) 
and {V\L) are both steeper than before. The exact val- 
ues of the change in slope depend upon our choices for 
the intrinsic model, for reasons we describe in the next 
subsection. 

The parameter choices given above, with the selection 
described above, result in zero-points, slopes (given at 
top of each panel) and scatter (0.22 dex, 0.33 dex and 
0.05 dex in the top left, right, and bottom panels, respec- 
tively) which are in excellent agreement with those in the 
Haring & Rix (2004) compilation (see Appendix A). Fig- 
ure m shows that residuals from these relations exhibit 
similar correlations to those seen in the real data (com- 
pare Figure IA21 residuals from the M, — a relation do 
not correlate with cr, L or L — a residuals, so we have not 
shown these correlations here). This agreement is non- 



trivial: for example, it is more difficult to match all three 
observed correlations if the selection is based on {V\L) 
rather than {L\V). Changes in the parameter values of 
more than five percent produce noticable differences. 

The correlations shown in the two right hand panels of 
Figure 2] are similar in the full and in the bias selected 
sample, suggesting that the selection has not altered the 
shape of this correlation. This suggests that the cor- 
relations seen in the right hand panel of Figure IA2I are 
real-they are not due to selection effects. We comment in 
the final section on what these correlations might imply. 

The result of using the biased {M,\L) relation (both 
slope and scatter) in the mock luminosity function to 
predict black hole abundances is shown as the dashed 
line in Figure [5l It overpredicts the true abundances in 
the mock catalog, which are shown by the solid line. The 
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Fig. 5. — Comparison of the intrinsic (solid) and selection-biased counts based on using velocity dispersion (dotted) and luminosity 
(dashed) as the indicator of M, . Hashed regions show the abundances infered by Tundo et al. (2006) by combining the M, — L and M, — u 
relations in the compilation of Haring & Rix (2004) with the SDSS-based luminosity and velocity functions (Bcrnardi et al. 2003b; Sheth 
et al. 2003). 



dotted line shows that the analogous cr— based predictor 
is, essentially unbiased. This is not surprising in view 
of the results shown in Figure [5] — the observed (M, |L) 
relation is biased towards predicting more massive black 
holes for a given luminosity, whereas the (M, |cr) is al- 
most completely unbiased. The hashed regions show the 
L— and cr— based predictions derived by using the rela- 
tions in the Appendix (taken from Tundo et al. 2006) 
with the Bernardi et al. (2003b) luminosity function and 
the Sheth et al. (2003) velocity function. They bracket 
the results from our mock catalog quite accurately. Note, 
however, that in our mock catalog, the true abundances 
are given by the solid line: the much larger abundances 
predicted from the luminosity-based approach are due to 
the selection effect. For reference, this intrinsic distribu- 
tion is well-described by equation ([H]). 

3.2. Toy model of the bias 

If current black hole samples only select objects which 
have logcr which is exactly n-standard deviations above 
the {V\L) relation, then 

V~{V) _L-{L) 



rvL + n\ I 



ay OL 

and inserting this in the expression above yields 



(5) 



(M.|L) (A/.) , T,v~r,i,ryi^ , L - {L) 

- -n , 1 r.L. (6) 



cr. 



VL 



CTL 



The slope of this relation is still a,|i; because the selec- 
tion was done at fixed L, the slope of the observed (M, |L) 
relation is not biased from its true value. However, the 



zero-point may be strongly affected. The magnitude of 
the offset depends on how much of the M, — a relation 
is due to the a — L and M, — L relations. If n is posi- 
tive (as Figure [T] suggests) and r,y > r,LrvL, then the 
observed relation predicts larger black hole masses for a 
given luminosity than would the true relation. The asso- 
ciated (M, |cr) relation is obtained by transforming L to 
V using equation ([5]): 



{M,\V) (M,) , V-(V) 



cr. 



cr. 



-nr,v 



■ r,v 

cry \r,vrvL 
' ir,L/r,vrvL) 



(7) 



' VL 



Both the slope and zero-point of this relation are affected, 
and the magnitude of the effect depends on whether or 
not r,L > T^vTvL- If the correlation between M. and 
L is stronger (weaker) than can be accounted for by the 
M, — V and V — L relations, then the observed {AI,\V) 
has a steeper (shallower) slope than the true relation, and 
a shifted zero-point, such that it predicts smaller (larger) 
masses for a given cr than the true relation. Note that 
this means it is possible for both the L— and cr— based ap- 
proaches to overestimate the true abundance of massive 
black holes. 

If the selection chooses objects with small L for their 
cr, then the resulting relations are given by exchang- 
ing Ls for Vs in the previous two expressions. In this 
case, {V\L) and (M,|L) both have biased slopes and zero 
points, whereas (Af. |cr) has the correct slope but a shifted 
zero-point. Once again, the sign and magnitude of the 
bias depends on how much of the M.— observable corre- 
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lation is due to the other two correlations. For the choice 
of parameters used to make the Figures in the previous 
subsection, r,vfvL — 0.69 « r.^: the M, — L relation 
is almost entirely a consequence of the other two cor- 
relations. Thus, this toy model shows why the (M, |(t) 
relation shown in Figure [3] was unbiased by the selection 
(which was done at fixed a). 

4. DISCUSSION 

Compared to the ENEAR and SDSS-B06 early-type 
galaxy samples, the a — L correlation in black hole sam- 
ples currently available is biased towards large a for a 
given L (Figures [lH2]) . If this is a selection effect rather 
than a physical one, then current determinations of the 
M, — L and M, — a relations are biased. We provided 
a simple toy model of the effect which shows how these 
biases depend on the intrinsic correlations and on the 
selection effect. 

To illustrate its use, we constructed simple models of 
the intrinsic correlation and of the selection effect. We 
then identified a particular set of parameters which we 
showed can reproduce all the observed trends: the corre- 
lations between M, and L or cr, the a — L correlation in 
black hole samples, as well as correlations between resid- 
uals from these fits (Figures [3] and [4]) are all reproduced. 
We do not claim to have found the unique solution to 
this problem. However, if we have correctly modeled 
the selection effect, then our analysis suggests that the 
bias in the (M, |cr) relation is likely to be small, whereas 
the {M,\L) relation is biased towards having a steeper 
slope (Figure [3|) . This has the important consequence 
that estimates of black hole abundances which combine 
the {M,\L) relation and its scatter with an external de- 
termination of the distribution of L will overpredict the 
true abundances of massive black holes. The analogous 
approach based on cr is much less biased (Figure [S]) . 

The choice of parameters which reproduces the 
selection-biased fits reported in Appendix A has intrinsic 
distribution 4'{M,) that is well described 



(/.(M.)dA'f. 



exp(-/)/3 



AM. 
M. 



(8) 



where 0» = 0.002/Mpc^ (a,/3) = (2.0,0.3) and = 
M,/M*, with M* a characteristic mass. This mass is 
related to the mean mass lO^M© by (M.) = M* T[{a + 
l)//3]/r(a//?), so M: = IO^Mq. Since this 0(Af.) is 
based on the SDSS-B06 luminosity and velocity func- 
tions, it is missing small bulges, so it is incomplete at 
small M,. 

The intrinsic correlations in the sample are 



M,\Mr ) = 9.16 - 0.41 {Mr + 22) 



(9) 



and 



M.I logo-) = 8.85 + 3.82 log(CT/200 km s"^) (10) 



with scatter of 0.38 and 0.23 dex, respectively. A bisector 
fit to the M, — L relation would have slope {cr,/(TL){r,L + 
1/^»l)/2 = —0.62; a similar fit to the M. — V relation 
would have slope {a,/av){r,v + ^/r,v)/'^ = 4.49. Notice 
that the slopes satisfy a.|L ~ a,\Y av\L- This is because, 
for this choice of parameters, the correlation between M, 
and L is almost entirely a consequence of the other two 
correlations. 

If we have correctly modelled the selection effect, and 
we do not claim to have done so, then our analysis sug- 
gests that the intrinsic joint distribution of M., lumi- 
nosity and a is similar to that between color, luminosity 
and cr: of the three pairwise correlations, the two corre- 
lations involving cr are fundamental, whereas the third 
is not (Bernardi et al. 2005). Since residuals from the 
(cr|_L) relation are correlated with age (Bernardi et al. 
2005), the correlations shown in the top right panels of 
Figures [4l and \X2\ suggest that, at fixed luminosity, older 
galaxies tend to host more massive black holes than less 
massive galaxies. Understanding why this is true may 
provide important insight into black hole formation and 
evolution. 
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TABLE Al 

Parameters of the black hole sample used in this work. HR=Haring & Rix (2004); MH = Marconi & Hunt (2003). 

is from observed absolute MAGNITUDE; Mr IS INFERRED FROM — Mbulgc RELATION. 



aHI^ Logio AC ^HR-M,„,,, ^MH ^^MH 

[km s~^] [Msun] [mag] [mag] [km s~^] [mag] 
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APPENDIX 

THE a-L CORRELATION 

The main text compares the a—L relation for the bulk of the early-type galaxy population with that for the black hole 
sample. This was done by using luminosities transformed to the SDSS r— band and scaled to = 70 km s^^ Mpc^^. 
Specifically, we convert from B, V, R, I and K-band luminosities to SDSS r- band using 

B-r=1.25, y-r = 0.34 r-R = 0.27, r- I =1.07, and r - K = 2.8; 

these come from Fukugita et al. (1995) for SOs and Es. We describe the samples and required transformations below. 

The black hole sample 

The main properties of the black hole sample shown in Figure 1 are listed in Table 1 . The sample was compiled as 
follows. We selected all objects classified as E or SO by Haring & Rix (2004) in their compilation of black holes. We 
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Fig. A1. — The M, — L and M, — a relations in the su bset of the Ha ring . 
hnes show the fits to the filled circles given in equations \A2\ and (IA3I I. 



Rix (2004) sample which wc describe in the main text. Solid 



excluded NGC 4342 as they recommend, as well as NGC 1023 which Ferrarese & Ford (2005) say should be excluded 
because it is not sufficiently well resolved. Haring & Rix provide estimates of the bulge luminosities and masses of 
these objects; these assume Hq = 70 km s"^ Mpc~^. For consistency we also rescaled the black hole masses from 
i?o = 80 km s~^ Mpc~^ to Hq = 70 km s~^ Mpc~^. We convert from the V, R and I-band luminosities the report 
to SDSS r using the transformations given above. We use the bulge masses to provide an alternative estimate of the 
bulge luminosities by setting 

login (Afbuigc/Mp) - 11-52. 



0.492 



(Al) 



Table 1 also lists estimates of the K-band near-infrared 



Mr = -22.34 

this was obtained from a linear regression of on Mbuigo- 
luminosity from Marconi & Hunt (2003) which are shown in Figure 2. 

A word on the velocity dispersions is necessary. Haring & Rix (2004) and Marconi & Hunt (2003) use the velocity 
dispersions provided by Kormendy & Gebhardt (2001) (two galaxies in Marconi & Hunt have a different value). The 
two estimates correspond to different apertures {r^ and central, which probably means re/8). Tremaine et al. (2002) 
provide a detailed discussion of aperture effects. However, it is not clear that their assumptions about the scale 
dependence of the stellar velocity dispersion are consistent with more recent measurements (Cappcllari et al. 2006). 
Since the difference in aperture is not sufficient to account for the discrepancy between black hole samples and the 
normal early- type galaxy samples we have not made any correction for this difference. 

The M, — L, M, — cr, and a — L correlations in the Haring & Rix (2004) sample, with luminosities estimated from 
Mbuige are 

'\o^M,\Mr\ = (8.57 ± 0.10) - (^■^Q -^ O-IO) + 22) (A2) 



with scatter of 0.33 dex, 

<^log M. I log cr^ = 

with intrinsic scatter of 0.22 dex, and 



.21 ±0.05) + (3.83 ±0.10) log, 

, 200 kms 



(2.41 ± 0.10) -(^:^1±«(M. + 22) 
2.5 



(A3) 



(A4) 



(Tundo et al. 2006). Figure [XT] shows the M, — L and A/, — a relations and these fits. The dotted line in Figure[l]of 
the main text shows equation (|A4|) . 

Figure IA2I shows that residuals from the M, — a relation do not correlate with residuals from the a — L correlation 
(top left), whereas residuals from the M, — L relation do (top right). Residuals from the Af, — a relation do not 
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Fig. A2. — Residuals from the M, — a relation do not correlate with residuals from the a — L correlation (top left); whereas residuals 
from the M, — L relation do (top right). Residuals from the A/, — L relation also correlate with residuals from the L — u relation (bottom 
right), but they correlate weakly if at all with a (bottom left). 



correlate with cr, L ov L — a residuals, so we have not shown these correlations here. However, residuals from the 
M, — L relation do correlate with residuals from the L — a relation (bottom right) and they may also correlate weakly 
with cr itself (bottom left). 

The top right panel suggests that, at fixed luminosity, objects with larger M, tend to have larger cr. The bottom 
right panel suggests that these objects also ten d to b e fai nt for their a. But whether or not these are entirely selection 
effects is unclear; this is because the relations (jA2| ~ (jA4|) are a combination of an intrinsic correlation and a selection 
bias, so that correlations between residuals are also combinations of true and selection-biased relations. 

Fits to early-type galaxy samples 

The velocity dispersions for the early-type galaxy population were scaled from their measured values to those they 
are expected to have within a circular aperture of radius re/8. For ENEAR, when the observed redshift is used as the 
distance indicator (i.e. no correction is made for peculiar velocities), 

(logfxlM,) = (2.300 ± 0.110) - ("-^^^ Q-"^^) (Mr + 22). (A5) 
\ / 2.5 
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Correcting for peculiar velocities leads to a biased relation, because the correction makes L depend on a. This biased 
relation is 



The galaxies in the SDSS-B06 are sufficiently distant that the SDSS-B06 a — L relation is not affected by peculiar 
velocities. It is 



Expressions (|A5|) — (jA7[) are shown in Figure [T] of the main text. They are taken from Bernardi (2006), where a 
discussion of the systematics in the magnitudes and velocity dispersion in the SDSS database can be found. 

Analogous quantities in the K-band 

The main text also studied the cr — Lk relation in the black hole compilation of Marconi & Hunt (2003). This 
relation is 



where the fit is made only to the objects which are also in the Haring & Rix (2004) compilation. Matching SDSS-B06 
to 2MASS, and so estimating a — Lk for early-types yields 




(A6) 




(A7) 
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